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1 Introduction 



There is a set of jet-shape observables that contain 1/Q power correction terms [1-3]. These 
non-perturbative contributions originate from soft-gluon emission with energies of the order of 
the confinement scale. 

The standard technology for determining these power terms relies on the introduction of a 
fake gluon mass. In a more systematic approach, it is the variable m^, entering in the dispersive 
representation for the running QCD coupling, that plays the role of the power correction trig- 
ger [4-6]. Within this approach the magnitude of the power contribution is expressed in terms 
of a standard dimensional integral of the non-perturbative part of the QCD coupling ("non- 
perturbative effective coupling modification"), multiplied by an observable-dependent (easily 
calculable) coefficient. 

In both cases the gluon decays, essential for building up the running coupling, were treated 
inclusively. This is known as the naive approximation. It is based on replacing the actual 
contribution of the final-state partons by that of their massive parent gluon. This procedure 
is intrinsically ambiguous since there is no unique prescription for including finite-m^ effects 
into the definition of the shape variable. Moreover, Nason and Seymour [7] questioned the 
application of the inclusive treatment of gluon decays. They pointed out that jet shapes are 
not inclusive observables, and therefore that the configuration of offspring partons in the gluon 
decay may affect the value of the power term at next-to- leading level in Og, which a priori is no 
longer a small parameter since the characteristic momentum scale is low. Both these problems 
can be resolved only at the two-loop level. 

For the thrust case the problem was addressed recently in [8]. The two-loop analysis led to 
a modification of the naive prediction for the coefficient of the 1/Q term by a perturbatively 
calculable Ca- and rij-dependcnt numerical factor, which has since been named "the Milan 
factor" [9]. A priori one might expect the Milan factors to differ for different observables. 

In this paper we consider, in addition to thrust, the following infrared- and collinear-safe 
observables: the C-parameter, invariant jet masses and jet broadening. In near-to-two-jet kine- 
matics (soft limit) these observables are linear in secondary particle momenta, and therefore 
produce 1/Q terms. 

We compute Milan factors for these observables, which rescale the naive predictions. We 
show that given an appropriate, consistent, definition of the naive approximation, the observ- 
ables belonging to the same 1/Q class possess a universal Milan factor, M. ~ 1.8. 

We address the problem of constructing renormalon-free answers that combine the trun- 
cated perturbative series with the power-behaving non-perturbative contributions. We choose 
to quantify the magnitudes of the non-perturbative contributions by means of finite-momentum 
integrals of the running QCD coupling rather than in terms of the moments of its disper- 
sive partner, the so-called "effective coupling" acfr [6], as suggested in [2]. The corresponding 
numerical "translation factor" 2/7r happens to practically compensate the effect of the Milan 
factor. As a result, we expect the phenomenology of 1/Q power contributions to jet shapes to 
be affected only slightly. 

The same Milan factor applies to the 1/Q power correction to the energy-energy correlation 
measure in e'^e~ annihilation (EEC) outside the region of back-to-back kinematics. A qualitative 
analysis of EEC in the special back-to-back kinematical region shows that one can expect large 
non-perturbative corrections characterised by fractional powers Q^^ , with 7 < 1. This results 
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from an interplay between non-perturbative and perturbative effects, which is typical for semi- 
inclusive observables. A fractional power should also be present in the distribution of heavy 
lepton pairs with small transverse momenta (the Drell-Yan process). 

The paper in organised as follows: in section 2 we define and describe the observables in the 
soft radiation limit; section 3 is devoted to the all-order resummation programme; in sections 
4 and 5 we extract the leading 1/Q power corrections to the radiator and to the corresponding 
distributions; section 6 is devoted to the merging of the perturbative and non-perturbative 
contributions in a renormalon-free manner; we summarise the results and present conclusions 
in section 7. 

2 Kinematics and observables 

Here we consider the kinematics of multi-particle ensembles consisting of a primary quark and 
antiquark with 4-momenta p and p, and ra secondary partons ki. Hereafter we shall consider 
the phase space region of soft secondary partons, soft gluons and their decay products, so that 
the primary p and p belong to the opposite hemispheres. 

We introduce the standard light-cone Sudakov decomposition in terms of two vectors chosen 
in the direction of the thrust axis: 



For a massless particle q, we have aqf5q = g|. In the soft approximation, for the primary quarks 
we have 



while all the secondary-parton momentum components are small. The differences 1 — ap and 
1 — [3p are linear in gluon momenta and will be taken care of; the components (5p and ap are 
quadratic and will be neglected in what follows. 

The jet shapes defined with the use of the thrust axis are T, the invariant squared mass of all 
the particles in one hemisphere ("jet mass" M^) and the sum of the moduli of the momentum 
components transverse to the thrust axis (the jet broadening B [10]). 

In terms of Sudakov variables they read: 



pi^ + pi^ = Qt'- = p'^ = 0; 2{PP) = Q^ = 1. 



The parton momenta are: 



p = apP + PpP + pt , 
p = apP + PpP + Pt ; 
ki = aiP + PiP + kti. 



(l-ap), (l-Pp) «1, 



(2.1) 



N 




(2.2) 



i=l 




(2.3) 



(2.4) 
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Here we have defined the right hemisphere as containing the quark p, so that for a particle in 
the right hemisphere, a < (3 (correspondingly, partons with j3 < a contribute to M£, Bl). We 
have also used the fact that the total transverse momentum of particles in each hemisphere is 
zero, which is a property of the thrust axis. 

There are two types of invariant mass and broadening distributions under consideration. 
One can discuss the total squared mass, the total broadening, etc., which are characteristics of 
the whole event: 

M| = M| + M| , Bt = Br + Bl. 
Alternatively, one can study the heavy-jet mass, the wide-jet broadening: 

= max {M^,M|}, = max{5R, S^} . 

The all-order perturbative analysis has been carried out both for the total and for single-jet 
("heavy", "wide") characteristics [11,12]. As we shall see later, the non-perturbativc effects in 
the single-jet characteristics, in general, amount to half of those in the total event observables. 
The C-parameter and the energy-energy correlation function (EEC) are defined as 

C = -j^PaPhSa? Oah , (2.5) 

ab 

EEC(x) = ^^a£;6(5(cosx-cose„b), (2-6) 

ab 

where Pa and are the modulus of the 3-momentum and the energy of parton a (the quark, 
the antiquark or one of the secondary partons); 9 ah is the angle between the two momenta. The 
sums over a, h run over all particles in an event. (Notice that within the definition of the C 
parameter and EEC, each pair a, h is counted twice.) In the massless particle approximation 
one has the sum rule 

C=\f^ d{cos x) EEC(x) sin^ x . (2.7) 
Now we consider each of these observables in the soft limit. 



Thrust. If only one gluon is emitted, say in the right hemisphere, ai < (5i < 1, then 

l — T = a\-\-ap-\-j3p. 
Since the thrust axis is such that = P^, we have 



p1 p1 a n 

Pi J- — Pi 

The differential thrust distribution takes the form [13] 

(l-T)^ = ^f41nJ--3 + 0(l-r) 

Here the logarithmic piece comes from the soft limit neglecting the quark contribution: ap = 
Pt ^ ai = Pt/Pij so that 1 — T ~ ai. The constant term accounts for gluon and quark 
contributions in the region of relatively hard gluons: /3i ~ /3p ~ 1, ai ~ ap. 
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In higher orders the soft gluon contributions exponentiate into the double-logarithmic Su- 
dakov form factor, while the effects of hard gluons and of quark recoil remain subleading, being 
down by one power of ln(l — T). However, the quark contributions to thrust, Up + (3p, don't 
contribute at all to the 1/Q power-suppressed contribution, which is driven exclusively by soft 
gluons. Therefore hereafter we ignore these quark contributions to (2.2), though we note that 
the perturbative part of the answer takes them into proper account [11]. 

Defining 

A = ^ai, A = ^Pi, 

and neglecting the quark contributions we have 

1-r = A + A. (2.8) 



Jet mass. It is common in the discussion of jet masses, in addition to considering the sum of 
the squared masses of the two hemispheres, 

Ml = M| + Ml , 

to study the heavy jet mass, 

M^ = max{M|,M|}. 

To simplify the analysis we postpone the consideration of which hemisphere is heavier and 
discuss the double-differential distribution in M| and M|. This will help avoid the considerable 
confusion that arises in a direct analysis of the non-perturbative effects in M|. Indeed, 
aiming at power effects, we should treat the two hemispheres on an equal footing, since it is the 
peHurhative radiation which determines which of the two jets is heavier. 

The quark contributions to the jet mass (2.3) can be analysed and dropped for the study of 
the 1/Q power corrections, as in the thrust case, so that we have 

M| = A{1 -A), Ml = A(l - A) . (2.9) 

In the first order in ag the relation between the jet masses and the thrust reads 

1 - T ~ M| M| = M| . (2.10) 

We observe that at the level this relation becomes non-linear. Prom (2.8) and (2.9) one 
derives 

,1_T,1±I = MJ.M2-M^. (2.11) 

This introduces corrections to (2.10) of relative order M| ^/Q^. In the region of small jet masses 
these effects are negligible. At the same time, averaging over, say M| in the M|-distribution 
will produce a cross-talk correction at a relative level of ^M^/Q^^ ~ asiQ"^)- As a result, the 
power term 1/Q in which we are interested will also acquire a relative correction of order as{Q'^). 
This subleading perturbative rescaling is beyond the accuracy that we consider. 
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Broadening. The distribution of the jet broadening, B, is different in two respects from 
those of the other variables considered here. At the perturbative level, the quark recoil can 
not be neglected. In the first order, for example, the quark contribution to B is simply equal 
to that of the gluon. Hence, the factor of two in the definition of B (2.4). In higher orders 
an account of quark recoil becomes more complicated [12, 14]. At the non-perturbative level, 
broadening has an extra coUinear enhancement, leading to a power correction proportional to 
{In Q)/Q. Broadening is the only observable currently under discussion with such a type of 
leading power correction. One can expect analogous (i.e. In Q-enhanccd) confinement effects in 
the Ef distribution in hadron-hadron collisions, in the oblateness measure, etc. 



C-parameter. In the soft limit, that is in the linear approximation in gluon energies, the 
relevant contributions to C come from counting (twice) the gluon-quark pairs. 



= 5.2 



Y,P°kf sin' Oi, + Y.f''? sin^ [ ■ (2-12) 
, i i ) 



The gluon-gluon contributions are quadratic in gluon energies and therefore negligible. The 
quark-antiquark contribution can also be neglected since sin^ O^p is proportional to the total 
gluon transverse momentum squared, and thus is also quadratic in gluon energy. Setting p° = 
p° = I and making use of the Sudakov variables, = \{ai + Pi) and sin^ 6ip = sin^ = 
4:aiPi/ [ui + A)2, we get 



Energy-Energy Correlation. The problem of power effects in energy-energy correlations 
has two aspects. In the region of finite angles, x ~ 1, the leading 1/Q power term originates 
from triggering a gluon with an energy of order of the confinement momentum scale. This 
contribution does not involve gluon resummation and can be obtained within a fixed order 
analysis. It is straightforward to derive the EEC function in the soft approximation: 

EEC(x) = -|^i:^-«*('-^)- p.") 

As we shall see later, the power correction explodes at sin x where multi-parton resumma- 
tion effects become essential. 



3 Observables and the resummation programme 

In what follows we concentrate on observables which require all-order soft gluon resummation 
to be carried out at the two-loop level. We introduce the label F as a general name for the 
observables 1 — T, B, C and the double jet-mass, so as to treat them simultaneously. For the 
broadening we consider the observable B = Bt- EEC away from the back-to-back region doesn't 
require resummation, and will be addressed at the end. 
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3.1 Observables 



We consider the F-distributions in the region F ^ 1. The observables 1— T, C and B are given 
as a sum of contributions v{ki) from each emitted parton ki. They can be expressed in terms 
of the multi-parton emission distribution dan by the following relations 

a{dlnV) dlnV ^ 

„ , / ra \ (3.1) 

m \ i=l J 

with m the total number of secondary partons in the final state (soft gluons or their decay 
products). We have V = 1 — T, C and 2B (notice that B is defined in (2.4) with a factor 2, as 

a half-sum of the moduli of the transverse momenta). 

For small V, which corresponds to all the final-state partons having small transverse mo- 
menta, one can write dam as a product of two factors, 

= C{as) dwn ■ 

a 

The first is an observable-dependent coefficient factor which is a function only of as{Q). The 
second is the observable-independent "evolutionary exponent" , which describes the production 
of n small-fct gluons off the primary quark-antiquark pair and their successive splitting into 
m> n final partons. Hence 

I{V) = C{a,{Q)) S(F, a,) , S(F, a,) = ^J dwn ® " £ ^('^^)) " (^"2) 

We note that dwn implicitly contains a sum over m in the range n < m < 2n. The soft 
distribution dwn is normalised to unity: J2nl '^'^n = 1, where the transverse momentum of 
each gluon is integrated up to Q. The independent emission distribution is valid only for small 
kti- It simplifies the treatment but mistreats the non- logarithmic region of large transverse 
momenta, kt ~ Q, both in real and virtual terms. This is compensated by the factor C{as) in 
(3.2). For small V values this factor is constant. The dependence on InV is embodied into the 
S factor. 

The essential momentum scales in the coupling in S range from VQ to Q. The perturbative 
treatment requires VQ S> Aqcd- The power correction we are interested in originates from 
v{ki)Q ^ Aqcd- 

The double-differential jet-mass distribution can be obtained via the distribution in A, A 
(see (2.9)) as 

^-'da _ f^^^-d^j:iA,A) 



dMj^dMl 



C{am) I dAdA ^-^^ ml - A{1 - A)) 5{mI - A{1 - A)) . (3.3) 
The double-distribution T,{A,A) is given by 

E(^, A) = Y,j dwn e{A - ai) QiA - ^ A) • (3-4) 

We denote it as the D-distribution (with D standing for double). 
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3.2 Soft multi-parton ensembles 

To obtain the multi-parton soft-emission formula at the two-loop level it suffices to take into 
account a single splitting of each of the gluons radiated off the primary qq pair. It has the form 

/ ^0[«s(O) + 4.x(fc|)]-4C^ / dr.(^„fc.)(g)'iM^}. (3.5) 

This expression has the following structure. The first term describes the real emission of a soft 
massless gluon (the term with an ill-defined q;s(0)) and the two-loop virtual contribution to it 
(the renormalisation-scheme-dependent function x(^t ))■ The last term on the first line describes 
the gluon splitting into a gg or qq pair. The corresponding matrix element is given in [8]. 
Finally, the exponential factor on the second line stands for the total virtual correction factor 
which ensures the normalisation of dwn- The expression (3.5) embodies the production of up to 
m = 2n secondary final partons. 
The two-parton phase space is 



,^ /, , N -rrdaicfkti dad^ktd^qt , , 

d£2{kxM) = IT— = -— z{l- z)dz, 3.6) 

-■--^ a,: TT a TT TT 



where a = ai + a2, kf is total transverse momentum of the parent gluon, qt is the relative 
"transverse angle" of the pair, and z is the fraction ai/a. We have 

ai = za, a2 = {1- z)a , kt = hi + kt2, Qt = ^ - ' (^•'^^ 

We introduce the rapidity and the mass of the parent gluon by 



= z{l - z)q^ • a = ^Jkf + m?e-'^ , /3 = ^Jkf + m?e'^ . (3.8) 



The right (left) hemisphere corresponds to r/ > (ry < 0). 
In these terms the two-parton phase space reads 

dT2{ki, k2) = dm^ dr] dz , (3.9) 

TT 2tT 

where (p is the angle between kt and qt. 
3.3 S— distribution 

Taking advantage of the factorised structure of the multi-gluon matrix element (3.5) we introduce 
the source function 

u{k) = u{a,l3) = exp{—i'v{k)) 
to write in terms of the Mellin integral transform 

i=l / i=l 
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Then the S-distribution takes the form, for F = 1— T, C and 2B, 

^^^^ = / 2^^ ^'"^ ^"""'"^ ' ^"""'"^ = E / n uih) . (3.10) 

For the D distribution in (3.4) we have to introduce two Melhn transforms: 

^A,A) = [ e'^^^ I e'^^^e-^M . (3.11) 

J 2m vi J 2m V2 

The radiator 7?. is a functional of the source u{k) and has the fohowing two-loop expression: 

+ ACf J dr2ih,k2)[^y ^M\h,k2) [l-u{ki)uik2)]. (3.12) 

For the specific observables under consideration the source functions are: 

T: u{k)=e-'"^e{P-a)+e-''^e{a-P), 
D : u{k) = e(/3 - a)e-''"^ + Q{a - p)e-''^^ , 

C: u{k)=e-^<^\ c{k)=6-^, ^^'^^^ 

p + a 

B : u{k) = e-^V^i^e*^^°'* . 

In the source for the distribution in total broadening, B = Bt, we have introduced (see [14]), 
in addition to the Mellin variable ly conjugate to the modulus of the transverse momentum, 
the two-dimensional impact parameter b conjugate to the vector transverse momentum, and we 
write bkt = 6\/a^cos$. 

In (3.13) we have expressed the sources as a function of the Sudakov components a and /?, 

u{k) = u{a,P) , 

where, for a massless parton, af3 = k^. 
3.4 The radiator 

The various terms in TZ are ill-defined and only the sum has physical meaning. In what follows 
we cast TZ as sum of three contributions which, for our observables, are collinear and infrared 
finite. To this end we introduce a source corresponding to the parent gluon momentum u{k\ + k2) 
and split the source in the two-parton contribution on the second line of (3.12) as 

1-U{ki)u{k2) = [l-u{ki+k2)] + [U{ki+k2)-u{ki)u{k2)]. (3.14) 

The expression in the first square brackets defines the "naive contribution" which treats the 
gluon decay inclusively. The second contribution to the radiator we refer to as the non-inclusive 
correction. The parent gluon momentum k\ + k2 is "massive" : 

u{ki + ^2) = u{a, P) , a = ai + a2 , P = Pi + P2 ceP = ^1 + m'^ ■ 
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We note that the use of the sum of 4- vectors of the offspring parton momenta in (3.14) is 
somewhat arbitrary. We could have used any combination of the momentum components which 
satisfies v = V1 + V2 in the limits of collinear and / or soft branching. The non-inclusive correction 
would then change correspondingly. In other words, both the naive coefficient for the power 
correction and the Milan factor depend on the definition of the "naive approximation," while 
the final answer remains unambiguous. Our prescription has the advantage of preserving those 
very simple numerical coefficients that are known from previous one-loop calculations. 

Noticing that the variables kfi, kt2, z and m? are invariant with respect to Lorentz boosts 
in the longitudinal direction, while a is not, we conclude that the invariant matrix element 
does not depend on a. This allows us to perform the a-integration explicitly and represent the 
non-inclusive correction in the factorised form: 

UnM = ACf J drr?dkl ^ dz (g) ' ^M^ , (3.15) 
with the non-inclusive "trigger function" 

= — [ U{ki + k2) - U{ki)u{k2) ] . (3.16) 

The m^-intcgral converges, in spite of the singular behaviour of the matrix element, oc 1/m^, 
because the trigger function 0,ni vanishes in the collinear limit, Vm?. 
The radiator then takes the form 



n[u] = nr., + 4Cf £ ^ + x(fc?)) n^ikh 



+ ACf C drr? dkl dz ^ ^'^S' ^o(fef + m^) , (3.17) 
Jo 27r 2! 

where we have introduced the "naive" trigger function 

^oik'^) = r —\l-u(a, P = k'^/a)] . (3.18) 

The inclusive integral of the two-parton matrix element on the second line of (3.17) is given, at 
fixed kt, by [8] 



/ 



r 

Jo 



where, for A;f <C 1, we have replaced the actual upper limit of the integral, m? < 1, by 
< 00 since the integral is convergent in the ultra-violet region. 

The m^-integral here is ill-defined since it diverges in the collinear two-parton limit — > 0. 

The quantity as(0) on the first line of (3.17) is also ill-defined. It combines, however, with the 
PQ-term in (3.19) to produce the running coupling at the scale kf in the one-gluon emission, TZq 
(see below). 
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The virtual correction x(A;^) also contains a collinear divergence. It can be written down in 
terms of the dispersive integral 

with S a scheme dependent number. In the physical scheme, in which the coupling is defined 
as the intensity of soft gluon radiation [15, 16], one has »S = 0. In the MS scheme, for example 
we have 

The collinear divergences in % at /v,^ and in the logarithmic term of (3.19) at — > 0, 
cancel for the collinear /infrared-safe observables we are dealing with. The remaining correction, 
the "inclusive correction" TZin, is finite. 

Thus the final result for the perturbative radiator TZ at the two-loop level can be cast as a 
sum of three terms 

TZ = TZq + TZin + TZni ■ 

This regularisation procedure was introduced and discussed in detail for the case of thrust in [8] . 

In the following we study the structure of these contributions in terms of the so-called 
effective coupling a^s which was introduced in [6]. It is related with the standard QCD running 
coupling by the dispersive integral 

- ^^dm\'^fil. (3.22) 



fc2 Jo (m2 + it2)2 

Its logarithmic derivative is the spectral density for as', perturbatively. 

Naive contribution. The first contribution is the so-called "naive" contribution which emer- 
ges from the inclusive treatment, that is when the contributions to V of gluon decay products 
are replaced by that of the parent "massive" gluon. It incorporates only the term proportional 
to Pq in (3.19) which, together with as(0) in the Born term of the radiator, builds up the running 
coupling in the effective one-gluon emission. We have 

Using (3.23) and integrating by parts we arrive at 

= — rf- aeff(m ) ^ ^oih +m) 

Cp dm? 2\ o / 2\ 

= / — o-Q!eff(m ) i2o(m ) . 
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Inclusive correction. The second contribution is the "inclusive correction" which takes into 
account the logarithmic piece in (3.19) dropped in the naive treatment. If massless and massive 
gluons contributed equally to V, this logarithmic term would cancel completely against the 
second loop virtual correction x to one-gluon emission. The mismatch is proportional to the 
difference of the massive and massless trigger functions. Using (3.23) and integrating by parts 
we obtain 



8CfCa f dm^ aeff(»T- ) d f dkj kf{ki + m^) ^ 
Po J 'm 47r d In 7 kf + 
= no{k^ + m^) - no{k^) . 



(3.26) 



Non-inclusive correction. Finally, the "non-inclusive correction" in (3.15) accounts for non- 
soft, non-collinear gluon decays into two partons ki, k2, which kinematical configurations are 
mistreated by the "naive" inclusive approach, as was pointed out by Nason and Seymour [7]. 
Using (3.23) and integrating by parts we obtain 

4Cf dm^ aeff(m^) d \ f^^ d<t> ^ f^^ 2 1 , .2 o 1 foo7^ 

which form is suited for extracting the power correction. 

In the following we discuss the three trigger functions, Oq, f^m and O^j, and then the 
corresponding power terms, STZq, 6Tlin, STZni- We first consider the "Linear" observables T,D 
and C, and then the special case of broadening, B, (the "Log-Linear" observable). 

3.4.1 The trigger functions for the Linear observables {V = 1—T,D,C) 
Naive contribution. The "trigger functions" ClQ{m?) for our observables are: 



/•mm 

r : no = 2 dr/ (1 - exp {-mue-''}) , 

Jo 

D : no= d?7(l -exp{-mz/ie-''}) + (l-exp{-mi/2e"''}) , (3-28) 

Jo 
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C : 0,0 = I d?7 ( 1 — exp \ —mv 

Jin 



im 



cosh 77 



Each rio vanishes oc Vrn? in the small-mass limit, thus ensuring the convergence of the m^- 
integration in (3.25,3.26). 

The main, perturbative, contributions to (3.25)-(3.27) come from the logarithmic integra- 
tion region 1/z/^ < < 1. The expression (3.25) for TZq reproduces then the usual double- 
logarithmic perturbative radiator, proportional to Qg In^ t^. The inclusive and non-inclusive 
contributions (3.26) and (3.27) produce two-loop subleading corrections of the order of aglnv. 
Actually, the standard resummation programme can only guarantee exponentiation of the terms 
that contain a" In"* u, with m > n. Therefore, the perturbative corrections to 7?.j„ and 7?.„j do 
not belong to the "radiator" and are taken care of by the exact two-loop calculation (the per- 
turbative matching procedure) [11]. 

The genuine non-perturbative contributions originate from the mass integration region <^ 
In this limit, we can expand the exponents in the trigger functions (3.28) to first order in 
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m to obtain 

J^o(m2; V) = 5no + (m^) ; Sn*^^^ = p(^) • m , 

where the p-parameters are given by a standard rapidity integral depending on the observable 
under consideration: 

piT) =2v dr]e-'> = 2u , p^^^ =u^ + u^^ = u I drj — — = Sniy . (3.29) 

Jo J-oo cosh 7/ 

Inclusive correction. In the small-m limit we have {V = 1 — T,D,C) 

S^Z^ = P^""^ ■ {^\/k! + m^ - fct) (3.30) 

with the p-factors as in (3.29). 

Non-inclusive correction. In the small-m limit the trigger function Q^i in (3.16) gets sim- 
plified: 

= P^^^ ■ {hi + h2 - ^jkl + m?) . (3.31) 
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{V) 
ni 



In the linear approximation, each source u{k) produces a transverse mass factor (-\/ + for a 
massive gluon, and ku for a massless parton), the rest being the universal observable-dependent 
factor p(^) defined in (3.29). The broadening case is analogous, as discussed below. 

3.4.2 The trigger functions for the Log-Linear observable {V = 2B) 

The broadening measure should be treated separately since the large-distance contribution to 
B is enhanced by a Inm^ factor as compared with the "Linear" T,D and C. 

Naive contribution. In the case of broadening the source is independent of r/ for fixed (and 
small) transverse momentum, or rather transverse mass. As a result, an additional log-factor 
originates from integrating over the rapidity of a parton: 

pin 

B : ^0 = [l - e-'''" Jo{mh)) j dr] . (3.32) 

In the small-m limit we obtain ^ 

S^i^^ = um-ln — :r . 

m^ 

This result, however, can be improved by taking into account the non-soft contribution to the 
gluon radiation probability. This can be done by simply using the exact quark quark -|- gluon 
splitting function instead of the soft da/ a spectrum in (3.5): 



/ dri = 2 — 
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With account of this correction, the small-m limit of the broadening trigger function becomes 

6n'>!^'^ = z/m • fin ^ - = i/m • In . (3.34) 

\ m"' 2 J 

Note that the 6-space resummation does not affect the leading power correction: the expansion 
of the exponent exp{ibkt) in (3.13) lacks a term linear in kt after integration over the azimuthal 
angle Thus, a small transverse momentum of a gluon, and not that of the primary quark, is 
relevant for the 1/Q power contribution to the broadening. 

Inclusive correction. Given the definition of the inclusive trigger function (3.30), we substi- 
tute the naive trigger function for the broadening, (3.34), to obtain 



g-3/2 g-3/2 

k^ + w? In - kt In „ 

kf + kf 



(3.35) 



Non-inclusive correction. From the definition of the non-inclusive trigger function, (3.16), 
for the broadening we derive 

<50if ) = V [kt, + - ^k^t+m^) In ■ (3.36) 

4 Power corrections to the radiators 

To extract the leading power contributions we replace the cS'ective coupling acfr(^^) in the 
integrals for TZq (3.25), TZin (3.26) and TZni (3.27) by its non-perturbative component ("effective 
coupling modification") (5aeff(m^) and employ the small-m^ limit of the corresponding trigger 
functions, 50. 

The answers will contain the non-analytic in m? moments [17] of the effective coupling 
modification, 

^2p,, = 2^ ^ ('^ f In-^ ^ 5«efr(m2) , (4.1) 

with ji an arbitrary momentum scale. We remind the reader that the moments analytic in m? 
vanish, that is for p integer and g = [6, 18]. 



4.1 Linear observables {V = 1-T, D, C) 

Naive contribution. We start from STZq, the power correction to the naive contribution: 

TT J Til 

where we have restored the dimensional factor 1/Q. Here Ai is the p = i, g = moment of the 
coupling modification (4.1): 

Ai = Ai o = Tr-f 5oies{w?) . (4.3) 
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Inclusive contribution. We have to consider the following integral in (3.26): 



d In m? \Jq + m? 

We have replaced here the actual upper limit of the fc^ integration, = 1, by oo since the 
correction term is proportional to In and does not produce the leading 1/Q power contri- 
bution. Introducing the dimensionless integration variable x = kt/m we obtain the "inclusive" 
power correction in the form 

with 

f X dx 1/2/1 2\\ n c\c\f\ ^ A. 



= / ln(a;^(l + x^)) = 3.299-^. (4.6) 

/3o h (l + a;2)(Vr+^ + x) ^ ^ >' /?o ^ ' 

Non-inclusive contribution. To calculate the l/Q power term in the non-inclusive contri- 
bution in (3.27) we have to consider the multiple integral of the two-parton decay probability 

"""s^irii'^^i"*''"^^"^ (.„+fe-;^)}. (4.7) 

The linear mass-dependence can be factored out, leaving the convergent dimensionless integral 
which was calculated in [8]. The result reads 

^^(T,D,C) ^ ^'ji{C,D,c) ^ ^ 2/?-i (-1.227CA + 0.365C^ - 0.052n/), (4.8) 

where the three terms originate respectively from the soft-gluon, hard-gluon and quark parts of 
the matrix element (see [8]). 

The final results for T, D, C. Finally, the leading power contributions to the radiators of 
the r, D and C distributions are: 

^7^(^'^'^) = = <57^f^'^)A4, (4.9) 

with the p factors given in (3.29). This answer differs from the naive predictions T^q^^ defined 
in (4.2) by the universal "Milan factor", 

M = l + Tin + Vni = 1 + /?o~' (1.575Ca - 0.104n/) 
= 1.490 (1.430) for n/ = 3 (0) . 

We combine the Milan factor with the Ai moment to define a small parameter 

V ^ . (4.11) 

The non-perturbative contribution to the radiator becomes 6TZ^^^ = p^^^V- 
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4.2 The case of Broadening 

The answer for broadening^ contains two moments: the In Q-enhanced term proportional to Ai 
and a new log-moment (4.1) with p = = 1: 

= ^^2^Jo ^ m\n-jda,s{m^) . (4.12) 

Naive contribution. The naive power contribution to the S-radiator reads, see (3.25) , (3.34) , 



(■B) _ 



(4.13) 



Since the value A'^ is not apriori known, we can eliminate it by choosing fj, which then becomes 
a new non-perturbative parameter. 

Inclusive contribution. For the 5-distribution (4.4) becomes 

The answer reads, see (3.35) 

where 

, 2Ca P xdx 



/3o 



foo r r , 1 Ca 

/ 77—^ ^ lii^^ - Vl+x2 ln(l+x2) In(x2(l+x2)) = -22.751-^ . (4.16) 



Non-inclusive contribution. The non-inclusive power correction to the broadening radiator 
is obtained by using (3.36), the small-m part of the trigger function, in (3.15): 

where r^j is given by the dimensionless integral similar to that for Vni in (4.8) but with an 
additional logarithmic factor: 



r' ■ 

m 



J 2Tr J J TT m kf + 



Its numerical value is 

(4.808Ca - 0.884CA + 0.116n/) , (4.19) 



r' ■ 

ni 



where the three terms, as before, correspond to the soft-gluon, hard-gluon and quark contribu- 
tions. 



^In this revised version we wish to point out that the analysis of the jet-broadenings presented hero is super- 
seded by that given in Eur. Phys. J. direct C5 (1999) 1. While the previous section has been updated to account 
for a mistake of a factor of two for the non-inclusive pieces in the original version of the paper, this section has 
not been updated (this would simply involve the inclusion of an extra factor of two in r'„i) since it would in any 
case not represent the full answer. 
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The final result for B. The full power term in the broadening radiator is obtained by 
combining the three contributions. One reconstructs the Milan factor and finds 



2M 



^i(ln^-e) -(Mi+^i) 



(4.20) 



where 



3 2 + r'„ + r' 



2M 



1.692 + 



0.0765 



CA/3-0.0392n/ 



1.930 (1.692) for n/ = 3 (0) . (4.21) 



Here we have introduced the combination of the moments, i^'^+^i, which will be of convenience 
later on, see section 6. 

Defining also the parameter V' , we can write 



The new unknown parameter, V' can be traded for the scale of the logarithm: 



(57^(^) = 2z/ • :P In 

The scale is observable-dependent and reads 

Qb 



Q_ 
Qb 



In 



(4.22) 



(4.23) 



(4.24) 



5 Power effects in the distributions 

The leading power contributions to the radiators are proportional to the corresponding Mellin 
variable u, see (4.9), (4.23). Inserting this into (3.10) one obtains a shift of the corresponding 
perturbative distribution [19] proportional to the V parameter (4.11). 
We arrive, for small V = l-T,D,C,B &i 



da 
df 
da 
dldl 



(1-T) 

{A, A) 



da 



dT 
da 



{1-T-2V) 



da 
dC 
da 

Ib 



dAdA 
da 
dC 
da 

1b 



PT 

{A-V,A-V) 



PT 



(C - ^-kV) 

■ PT 

{B-VHQ/Qb)) 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



FT 



These expressions are valid as long as V , though numerically small, stays at the same time 
much larger than Aqcd/Q, that '\s V ^ V. The kinematical region V < V '\s dominated by 
confinement physics: the jets get squeezed down to small-mass (exclusive) hadron systems. 
By expanding the general answer for the distributions in the Taylor series in In V, 



f{V) = fMV-cvV) = fpT{V) - ^-^fUv) + \ 



(/pT ~ /pt) + 



(5.5) 
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with g'{V) = dg/dlnV, we observe that the true expansion parameter is V/V, rather than 
V. Such an enhancement at small is a common feature of the power contributions. For 

example, the l/Q^ power terms in DIS structure functions and e+e~ fragmentation functions 
are relatively enhanced near the phase space boundary as (1 — x)^^, where x is a usual Bjorken 
(DIS) or Feynman (e~^e~) variable [6]. 

Strictly speaking, within our approach only the leading power correction was kept under 
control. Therefore we should have dropped the rest of the series, starting from the quadratic term 
O {{V /V^Y spite of this we prefer to present the result in terms of the shifted perturbative 
spectra, because this form is easier to implement practically. The subleading power contributions 
become comparable with the leading one ioi V >V and remain to be studied. 



5.1 Jet mass distribution 

The invariant jet mass distribution (s) can be derived from the double differential distribution 
(5.2) with use of the relation (3.4). Generally speaking, the relation between the Sudakov 
variables A, A and the jet masses, M|, is non-linear: 



dAdA 



dMl dMl 
J 



J =1- A- A = ^1 - 2(M2 + Ml) + (Ml - M2)2 = l + O (M^) 

The relations inverse to (2.9) are 

A = Ml{l + Ml-e) A = Ml{l + Ml-e) 
^ = l + O (m2) . 



e = 



Ml 



Mf + J 



(5.6) 
(5.7) 

(5.8) 
(5.9) 



Calculating the mass distribution in terms of the shifted A, A perturbative spectrum in (3.4) 
amounts to shifting the values according to 



dMl dMl 



(Ml Ml) 



(fa 



dMl dMl 



{Ml 



VJ + V^, Ml-VJ + V"^) 



J PT 



The quadratic terms should be dropped within our accuracy. As for the linear power shift, in 
the small-mass limit, M|,M| <^ 1, we can approximate VJ ~ V. Therefore, in the small M^ 
region we can simplify the answer as 



^ ^(Ml,Mh 

dMidMr ^ 



d^a 



dMl dMl 



:{Ml-V, Ml-V) 



(5.10) 



PT 



Notice, however, that in the calculation of the mean squared invariant mass(es) the Jacobian 
factor will induce a perturbative correction to the non-perturbative shift of the order of V-as{Q^), 
since ( J) ~ 1 + (M^) = l + (as)- 



5.2 Power effects in whole-event and single-jet shapes 

Some shape observables characterise a single jet rather than the whole (2-jet) ensemble. Among 
these are such characteristics as the heavy-jet mass, the wide-jet broadening, etc., to be com- 
pared with the mean jet mass, and the total B. Within the first order analysis aiming at power 
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effects, tfiat very gluon tfiat induced tfie power-befiaving sfiift in tfie jet distribution, was mak- 
ing tfie jet under consideration "fieavier" ("wider") than its partner jet devoid of any radiation. 
This resulted in a strange unphysical picture in which all the non-pcrturbative effects were con- 
tained in the heavy/wide jet alone [1]. A proper treatment of the heavy-light relationship has 
been given by Akhoury and Zakharov in [20]. 

To illustrate the solution of the puzzle, let us take the total and heavy-jet mass distributions 
as an example. For the distribution in Mj, = -|- M| we immediately obtain the doubled shift 
of the corresponding perturbative expression, as compared with (5.10): 



— (M|) 



dMr 



) /"OO 

dMl / dMl 6{M^ - M| - Ml) 
Jo 



dM]^dMl 



{Ml -V,Ml- V) 



J PT 



-II 



dMl dMl 6{mI -2V - Ml - Ml) 



d^a 



dMl dMl 



(Ml, Ml: 



PT 



da 
dMl 



{Ml - 2V) 



J PT 



(5.11) 



To calculate the heavy-jet mass distribution one needs to perform the integration over the smaller 
of the two jet masses (for definiteness, we take Mh = Mr > Ml and double the answer): 



dMl 



^(Mi) = 2 j 

r 



dMl 



d^a 



dMl dM? 



{Ml -V,Ml- V) 



M%-P 

2 / dMl 



d^a 



PT 



dMl dMl 



{Ml - r, Ml) 



(5.12) 



PT 



da 

dMl 



{Ml-V) 



PT 



Recall that the perturbative answer for the integrated heavy-jet mass spectrum is essentially 
the squared single-jet distribution [11], 

Sj^(Mi)= (So„e-jet(Mi))', 

where Sone-jet describes a single-inclusive distribution in, say, M^, integrated over M| without 
constraint. This corresponds to putting 1/2 = in (3.11) and (3.13). Thus we conclude that while 
the total-mass distribution acquires the doubled shift, 2-V, {1^1 = 1^2 in the Mellin transformation 
language), the heavy- mass spectrum is sensitive to confinement effects in a single jet only {1^2 = 
0), and therefore the corresponding perturbative expression gets shifted hy 1 -V. 
The same pattern holds for the total- versus wide-jet broadening: 



da 

dBy/ 



(B 



w 



da 
dB]^ 



{Bw-'^-VHQ/Qb)) 



(5.13) 



PT 



to be compared with the twice larger shift in the total broadening spectrum (5.4). 



5.3 Energy-Energy Correlation 

At the perturbative level, the gluon-quark and quark-quark contributions to the energy-energy 
correlation function EEC(x) defined in (2.6) are comparable. As far as the power effects are 
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concerned, the leading 1/Q contribution to EEC originates from the correlation between one 
of the primary quarks and a very soft gluon, EEC^^. The gluon has an energy of order of the 
confinement scale, ~ /c^ ~ Aqcd- It is easy to see that this power correction to EEC acquires 
the same universal Milan factor as the T,D,C and B distributions: 



EEC(x) = EECPT(x) + ^ ^ = EECPT(x) + ^ • (5.14) 

sm'^ X Q siir x 



While we will not present here explicit results for the quark-quark contribution, EEC^^^, it is 
straightforward to show that it is of the form 

EEC^^P - ^2,1 



Both the quark-gluon and quark-quark contributions are more singular at small sinx than the 
perturbative part of EEC: 

EEC^/(X) A, EEC^P(x) ^2,1 

EEC^^ix) ^ Qsinx' EECPT(;^) " Q^sm^x' 

This, as we have discussed above, is a standard feature of power-suppressed contributions. The 
power-suppressed contribution from the gluon-gluon correlation may contribute as much as the 
quaxk-quark component, but it requires more study. 

In the region of very small sinx values, all-order effects in the perturbative quark-quark 
correlation become essential, which makes the perturbative correlation finite in the back-to- 
back configuration, that is for x = ti" [21,22]. Qualitatively, what happens is that the quark and 
antiquark are not exactly back-to-back, with the typical acollinearity angle, X) being given by 

sinx = (Q6)-^ (5.16) 

Here b = b{Q) is the characteristic impact parameter which determines the perturbative double- 
logarithmic distribution. A naive estimate gives: 

m - oc for nf = 3. (5.17) 

The acollinearity of the quarks leads to the freezing of the singularities in (5.15) when tt — x ^ X- 
For the back-to-back correlation we therefore expect 



EEC^/(vr) EEC^ 
EECPT(7r) ^^^^ ' l^C^ 



AibiQ), ^^^PT - A2,ib'iQ). (5.18) 



We expect the same fractional powers of Q for the non-perturbative corrections to the value 
of the energy- weighted particle flow in the current fragmentation region in DIS, in the photon 
direction (qt = 0). In addition, for the value of da/dM'^ dq^, the height of the plateau in the 
transverse momentum distribution of Drell-Yan pairs (W, Z bosons) at = 0, we expect a 
contribution which is quadratic in b{Q), 

exp (1^2,1 •6^(Q)}. 

PT 

A more detailed analysis reveals that the estimates for the values of the exponents are too naive, 
but that non-integer power effects are present [23]. 



da 




da 




dQ^dql 


qt=0 


dQ^dq^ 


qt=0 
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6 Merging the PT and NP contributions 

We now discuss how to obtain the fuh expressions for the distributions by merging the per- 
turbative and the non-perturbative contributions. We follow the procedure suggested in [2]. 
The answer we shall arrive at will differ, however, from that of [2] by a factor normalising the 
power terms. Part of it is the Milan factor we have discussed above. As we shall shortly see, 
an additional factor, 2/7r, arises due to the "translation" from the effective coupling acfj (which 
we use to trigger the large-distance power contributions within the dispersive approach of [6]) 
to the standard QCD coupling ag- 

For the radiators of the V = 1 — T,D,C distributions which contain 1/Q power terms we 
have obtained 

7^(Q2)=^PT(g2)^^(y)p^ (g_l^ 

where V is given in (4.11) in terms of the moment Ai of SacS- 

This representation is symbolic. Indeed, its perturbative part is given by an expansion which 
is factorially divergent. Meanwhile, the non-perturbative "effective coupling modification", da^s, 
also implicitly contains an ill-defined all-order subtraction of the pure perturbative part off the 
full effective coupling a^s- 

These two problems are of a similar nature. The uncertainty can be resolved at the price of 
introducing an "infrared matching scale", /ij. To this end we represent the answer as a sum of 
three terms: 

7^(Q) = 7^PT(Q;iV) +7^NP(Q,^J) - n^-^^%Q,fii;N) + O (af +1) , (6.2) 

with N the order at which the perturbative expansion is truncated. 

The factorially growing PT-expansion terms (infrared renormalons) in the first (PT) and 
the last ("merging") pieces in (6.2) cancel. The ///-dependence in the second and the third 
contributions approximately cancels as well, the cancellation improving for larger N. 

The structure of (6.2) is general; the same approach can be applied to incorporate the 
leading power contributions into different observables p, q. We proceed with the radiators for 
V = 1 — T,D,C distributions (p = = 0) taken as an example of the general matching 
procedure (6.2). To this end we split the running coupling, formally, into two pieces, 

a,{e) = a^^ik') + afP(A;2) . 

Here a^^(A;^) should be understood not as a mere function but rather as an operational pro- 
cedure: its momentum integrals should be evaluated perturbatively, that is order by order in 
perturbation theory as a series in as = as{Q^). Its non-perturbative counterpart, a^^(fc^), 
corresponds to the Saesim^) part of the effective coupling. It is supposed to be a rapidly falling 
function in the ultra-violet region [6, 24] . In the infrared momentum region neither of the two is 
a well-defined object. However, the full physical coupling as{k'^) is supposed to be finite down 
to k'^ ^0. Its infrared momentum integrals will play the role of phenomenological parameters. 
To define these parameters we first use the formal dispersive relation (3.22) for and its 
partner, Sa^s, to deduce 

/ dka^^{k^) = - ^mSaesim''). (6.3) 
Jo 4 Jq 



20 



Then, for Ai we have 

A^ = ^J^^ dk ar(fe') = ^ I" dk aTik") + O (^as'^^'Uf )) • (6-4) 

Here we have introduced a scale jii above which the coupling is well matched by its logarithmic 
perturbative expression, so that the error proportional to /Lt/a^^(/i|) induced by truncation can 
be neglected. Now we substitute 

into (6.4) and approximate Ai as 

Ai^^ r dk a.ik') - ^ r dk arik") . (6.5) 

Jo TT Jq 

The first contribution can be expressed in terms of a (/^/-dependent) phenomenological param- 
eter ao, 

/ dk as{k'^) = m ■ aoiiii) . (6.6) 
Jo 

Thus, the power contribution in (6.1) splits into two terms: 

The first term is the integral over the infrared region of the full coupling, 



(6.8) 



The subtraction term is given by the integral of the perturbative coupling, 

7^--s^(Q,A^,) = pW^:^ r dkoFJik^). (6.9) 



This perturbative subtraction contribution is obtained by substituting for (fe ) its pertur- 
bative series, 

00 

o^J^k^) = J2<Pe (in f) , as ^ as(A^I) . (6.10) 

Here fi^ is the renormalisation scale which hereafter we set equal to Q. The representation 
(6.7) is still symbolic, since the all-order perturbative subtraction diverges. Indeed, integrating 
(6.10) term by term one finds 



Jo 



00 ^ 

dk a^^ik^) = ^i-a^Y, (/?o^) Q , (6.11) 



i=o 



with Q being factorially growing coefficients. This behaviour is general. Taking, for the sake 
of illustration, the one-loop coupling, in which case (6.10) is a simple geometric series, 

«s ) = «s2^ ( '^0 2:^ln- 1 , A = —3 3^, (6.12) 
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one obtains 

= L T' dk( In^^ =f dte-H\ t = ln^. (6.13) 

Q Jo \ kj J\n{Q/ni) k 

For finite i we have Q ~ (ln(Q////))^, while for £ 3> ln(Q/ the coefficients start to grow as 
Ci ~ i\ (infrared renormalon) . 

The perturbative part of the answer, TZ^^ in (6.1), bears the same (infrared renormalon) 
divergence. As we have anticipated above, the factorial behaviours of the coefficients of these 
two perturbative expansions cancel. Indeed the corresponding coefficients in TZ^"^ for large i are 
given by integrals identical to (6.13) but with the upper limit taken as k<Q (t>0). Therefore 
the factorial divergence in TZ^"^ originating from the low-momentum integration region in the 
Feynman diagrams, k Qe^^', is subtracted off exactly by 7^™<'i'ge_ The mismatch vanishes 
at large £ S> hi^. Due to this cancellation, the difference {TZ^'^ — TZ'^^^^^) is well defined 
perturbatively, so we can truncate the perturbative expansions for both TZ^"^ and af^(fc^) in 

^^rnerge ^^^^ g^^j^g jyth ^^^^^ 

The merging piece truncated at £ = A'" — 1 takes the form 



N-l 



n--^%Q,„:-N) = p(^) ^ g as E (/^of^) Q • (6.14) 
The final answer reads 



n{Q) = 7^PT(g; N) + p(^) ^M^{ao-a,Y,{(^o^ycA+ ^niQ. f^i;N) . (6.15) 



N-l 



Q 



The answer contains two components, a perturbative piece with a smooth logarithmic Q- 
dependence, TZ^^ = TZ^'^ {as{Q^)) , and a steep power term TZnp — fZ^^^^^ oc m/Q. Corre- 
spondingly, the error A-ji of the representation (6.15) contains two separate components. 
For a moderate N the error in the perturbative part is obviously given by 

= 0{a^-^\Q')) . (6.16) 

The power-behaving component contains several sources. The first is the truncation error in 

•^merge ^l^jch amOUUtS tO 



^merge ^ ^ (^a^iQ^) = TZ''^ ■ O 



^A-), A./3of^ln^<l. (6.17) 

"0 / 27r /x/ 



The A parameter is smaller than unity provided that fij > Aqcd- It is important to stress 
that these estimates are uniform in A', since the factorially growing coefficients in TZ^'^ and 
•^merge cancel. Constructing the difference of the corresponding renormalon contributions, for 
the remainder of the series for TZ^^{Q, N) - 7^™^^s^(Q, iJ.i;N)we have 



The renormalon leftover being smaller than (6.17) prevents the magnitudes of the errors in 
(6.16) and (6.17) from exploding at large A. 
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Another error originates from the cutoff on the momentum integral in (6.4) at the finite 
value Hi: 



A^P = O (ar(/^l)|) = n^^-0 (^^) ■ (6.18) 

Within the logic of the present approach, is a rapidly falling function, so that by choosing 
a sufficiently large one can make this contribution arbitrarily small. 

The largest contribution to the error in the power-behaved piece will of course come from 
higher corrections to the Milan factor which will have an effect of the form 



where enters at some small scale. To quantify such effects one has to go beyond a two-loop 
calculation. Whether they are reasonably small depends on the effective interaction strength 
at small scales. The present-day phenomenological estimate of as/ it ~ 0.2 may lead one to be 
hopeful. 

To conclude our discussion of the errors, we note that the difference {JlF'^ — TH^^"^^^), the 
"regularised" perturbativc contribution to the observable, essentially corresponds to the in- 
troduction of an infrared cutoff in the frequencies in Feynman diagrams. The regularised 
perturbativc contribution, as expected, depends on the infrared cutoff jii. 

The /^/-dependences in TZ^^dij) and TZ^^'^^^^Q, m; N) compensate each other. With increas- 
ing N, the residual jLt/-dependence disappears. One can verify this by evaluating the derivative 
of (6.15) with respect to fij: 



din III 



n{Q) = O (A^^'-s" + A^P) . 



A final remark concerns the scheme dependence. Explicit expressions for the scries for both 
n^^iQ, N) and 7^™'=^g°(Q, fir, N) depend on the scheme, that is on the choice of the pcrturbative 
expansion parameter ag. Throughout our analysis we have been using the so-called physical 
(CMW) scheme [15,16] in which the intensity of soft gluon radiation equals, by definition, ag. 
Shifting, for example, to the popular MS scheme amounts to substituting 



a2_ 

«s ^ «MS+^^ + --- (6-19) 

in the pcrturbative series, with K as defined in (3.21). As for the non-perturbative parameter ao 
(and Oq, see below) it is defined in (6.6) in terms of the CMW coupling. Changing the scheme 
would change the expression in the left-hand side of (6.6), but would not affect the value of ao- 
Notice that the normalisation of the second, power behaving, term in (6.15) namely, 

7^^P(Q,/X7)-7^--g^(Q,Ai7;iV), 

differs by the factor 2M/Tr from that of [2] which is being currently used in the experimental 
analyses. We return to this point in the conclusions. 
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6.1 N = 2 merging 

For the practically important case N = 2, when the perturbative answer is known at the second 
order in the coupling constant as = as((5^), the relevant moments determining the power 
contributions to T,D,C and B radiators read 



2Cf 



27r 



In 



1 



and 



1 A' 

2 



(6.20) 



(6.21) 



Here the non-perturbative parameter ao has been defined in (6.6), and g'q is given by a similar 
integral but with an extra logarithmic factor: 

ao(Ai/) = / — "s(k), aQ[i2i) = / — as(fc ) In — . (6.22) 
Jo Jo 1^1 f^i 

Let us note that from (6.22) we expect the new non-perturbative parameter to be negative. 

Using the dispersive representation (3.22), it is straightforward to derive the following relation 



for + Ai 



1 a; + A,= '^ Tdk «^P(/c2) In - . 



(6.23) 



In analogy with the treatment of the Ai moment, (6.4), we truncate the integral on the right- 
hand side at the matching scale, fij: 

^ r dka^'^'ik') In- r dka^'^ik^) In-. 



Linear observables. Thus, the merging contribution to V = 1 — T,D,C distributions is, to 
second order. 



7^^;J^^(Q,^,;2)=p(^) 
and the final answer takes the form 
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7^ = 7^PT(Q;2) + pW.l^ 
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(6.24) 



(6.25) 



Log-Linear observable. Substituting (6.20) into the non-perturbative broadening radiator 
(4.20), we arrive at 
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merge 



2v 
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Q L 



ao ( In — - ^ 

Q 



a. 



27r 



Q 



as ln^-e+1 +/5o:t^ In' — + (2-0 In— + 1 
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The final answer takes the form 



7^ 



{B) 



n^^ + 2.^M^' 
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In 



1^1 



M (ao-as- Po 
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a[, + as + /3o^(ln-^ + 2 



27r 
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(6.26) 



(6.27) 



(6.28) 
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7 Summary and discussion 



In this paper we have considered the group of event shape observables that exhibit 1/Q non- 
perturbative corrections. We have shown that the two-loop effects result in a rescaling of the 
"naive" perturbative estimate of the magnitude of the power terms by a universal Milan factor. 
This is true for the non-perturbative effects in the thrust, in the invariant jet mass and in 
the C-parameter distributions, and for the main In Q-enhanced power term in jet broadening. 
The same Milan factor also applies to the energy-energy correlation measure away from the 
back-to-back region. 

The universality of the Milan factor for the linear jet shape observables in DIS has recently 
been demonstrated by M. Dasgupta and B.R. Webber in [25]. 

We have shown, following [19] that the 1/Q power effects in jet shape distributions result in 
a shift of the perturbative spectra, 



— {V) 



dV ^ 



cvV) 



V 



Ml 



PT 



Ml 



Q2' g2 



du 
dV 



(F) = 



da 
dV 



V - cvV In 



Q_ 

Qb 



V = Bt, Bw ■ 



(7.1) 
(7.2) 



PT 



The all-order resummed expressions for the perturbative spectra in (7.1) can be found, for the 
thrust and jet-masses, in [11], and for the C-parameter in [26]. For the perturbative broad- 
ening spectrum (7.2) one should use the recent theoretical prediction [14], which improves the 
treatment of subleading effects compared to the original derivation [12]. 

The answers are valid in the region 1 S> F S> 7^ ~ Aqcd/Q- The relative coefficients cy are 
given in the table 



V = 


1-T 
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M| 


Ml 


Bt 


Bw 


cv = 


2 


37r 


2 


1 


1 


1 

2 



(7.3) 



Note that the power shifts in the single-jet characteristics (M^, Bw) amount to half of those in 
the total distributions {M^,Bt)- 

The non-perturbative parameter V was defined in (4.11). The i?-distribution contains an 
additional non-perturbative parameter, the scale of the logarithm Qb, related to V' , (4.24), the 
first log-moment of the "coupling modification" SaeS- However, as was discussed in section 6, 
the parameters V and T" (or, Ai, A[) suffer from the same infrared renormalon ambiguity as 
the perturbative series. To remove this ambiguity one has to introduce the "infrared matching" 
momentum scale ///. At the two- loop level, the expressions for the shifts in terms of the non- 
perturbative ///-dependent phenomenological parameters uq (6.6) and Qq (6.22) are: 

>,2 



2 Q 



ACf 
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aoilJ'i) -as-/9op- 

ZTT 



2 Q 



2tt 



(7.4) 
(7.5) 



Here we have given the MS expressions, with K defined in (3.21) 
the shift in the broadening distribution (7.2) as 



Pin — 

Qb 



V 



The parameter V' enters into 



(7.6) 
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with the numerical parameter ^ defined in (4.21). Recall that we expect the parameter to be 
negative if the QCD coupling is to remain positive in the infrared region, see (6.22). The InQ— 
enhancement of the 1/Q contribution to jet broadening remains to be verified experimentally. 

By now there is a good experimental evidence for 1/Q power terms in jet shape variables 
both in e'^e~ [27] and DIS [28]. The magnitudes are consistent, within a 20% margin, for 
different observables. An estimate of the basic Ai parameter can be obtained, for example, 
from the second order fit to the Q-dependence of the mean thrust. 



with ai,a2 the known perturbative coefficients. If one boldly used the original relation (4.11), 
V = const/Q, the fit to the data would produce Ai ~ 0.15 GeV. However, such an estimate is 
meaningless: the very representation (7.7) would be renormalon-infested. Therefore one should 
use instead in (7.7) the representation (7.4) for V which is constructed as the difference of the 
pure power contribution, TZ^^ , and the perturbative subtraction, 7^™erge Thus, the price one 
pays for avoiding the rcnormalon ambiguity is the introduction of the finite infrared "matching 
scale" /^z/Aqcd which enters into the shifts (7.4), (7.5) both explicitly and via and Og. 
As we have argued above in section 6, the ///-dependence gets weaker when higher orders of 
the perturbative expansion are included. The accuracy of the merging of the perturbative and 
non-perturbative contributions is of order 



with the order of the highest perturbative term included (N = 2 in (7.7)). A reasonable value 
to take for /x/ is /x/ = 2 GeV. 

Having considered the problem of the infrared renormalon divergence of the perturbative se- 
ries, we have chosen to quantify the non-perturbative effects in terms of the momentum integrals 

of the running coupling as{k'^) rather than in terms of the moments of the non-perturbative 
effective coupling modification, dacs^m?), as was suggested in [2]. Our non-perturbative param- 
eter ao differs by a factor 2/7r from the analogous parameter ao of [2]. This factor originates 
from the "translation" relation 



The power contributions {1Z^^ —TZ^^'^^^) of [2,19] were written in terms of ao and did not include 
the two-loop effects (the Milan factor). Thus, to switch to the expressions of this paper (based 
on ao, a'o), they should be multiplied by 2/7r • A^. These two factors practically compensate 
each other: 2A4/ it ~ 1.14 (nj = 3). Strictly speaking, the remaining small renormalisation could 
affect the phenomenological fits to the data, since it multiplies both the non-perturbative term 
TZ^^ and the perturbative subtraction contribution (T^mcigc-)^ ^j^g latter possessing a residual 
log Q-dependence. 

We could have chosen to express the answer in terms of the moments of as suggested 
in [2]. This would have required, for the sake of consistency, the reformulation also of the 
standard pt series and the pt subtraction in terms of acfj. Carried out to all orders, this would 
have given the same answer. However, if we restrict ourselves to A = 2 we observe the following 
puzzle: neither TZ^^ nor 7^™<^>'g'^ is sensitive to the difference between as and a^s to this order. 



(1 - T) = aias + a2as + 2P 



(7.7) 
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As a result the coefficients of Tl'^^'^^^{aes) and Tl™^'^^^{as) differ by the factor 2/tt, apparently 
for no good reason. This mismatch can be absorbed, but only partially, if one chooses an 
appropriate relation between the cutoff scales /x/ in the two representations. We believe that 
the origin of the puzzle lies in the fact that the relation between as and aeff intrinsically contains 
the non-perturbative region, and so only holds when one considers all orders of the perturbation 
series. We intend to consider this problem in more detail in a future publication. 

Concluding the discussion of the Milan factor we should stress that only after having gone 
through the ordeal of the two-loop analysis, can one unambiguously predict the relative magni- 
tudes of the power-behaving confinement contributions to different observables. This is because 
the corrections to the magnitude of the power term which are quadratic in small-scale as, get 
promoted to a finite renormalisation factor. This can be understood with the help of the fol- 
lowing qualitative argument. At the perturbative level, the correction can be recast in terms 
of a rescaling of Aqcd in the argument of as in the one-loop answer. At the level of the power 
terms, however, such a rescaling of the argument of the coupling immediately triggers a cor- 
responding rescaling of the dimensionful numerator of the 1/Q contribution. In spite of the 
entirely illustrative nature of the argument, it gives a hint as to why we do not expect the 
higher loop corrections to be as crucially important as the two-loop one. Indeed corrections 
beyond the second loop do not affect the characteristic momentum scale of the problem and 
thus the normalisation of the power contributions acquires only a factor of the form 1 + (as), 
which is formally subleading. Whether the relative size of such higher order corrections can be 
kept under control depends on the actual strength of the effective interaction at small scales 
and remains to be quantified. 

Given the importance of the two-loop effects, the "naive" first order approach becomes 
intrinsically ambiguous. Within this treatment one introduces into the Feynman diagrams a 
coupling that runs with the gluon virtuality (fake gluon mass) in order to define a "projection" 
onto the confinement interaction region. Such a procedure is non-unique with respect to the 
precise determination of the argument of the coupling, and the inclusion of the finite gluon 
mass effects into the kinematics and into the definition of the observables. For example, in the 
definition of T for the {qq+ massive gluon) system suggested by Beneke and Braun in [18] , the 

— * 

finite mass effects were included in the normalisation of thrust {pq + Pq + \kg\ < Q) . Being as 
good a naive approximation as any other, this produces an answer which differs by a numerical 
factor from the "standard" one. Our claim is that if the two-loop effects were evaluated, the 
final result based on the Beneke-Braun prescription would coincide with the one obtained in [8] 
and generalised in the present paper. 

Beyond the naive treatment, one can no longer embody all gluon decays into the running of 
the coupling. Those decay configurations that affect the value of the observable under consider- 
ation cannot be treated inclusively [7]. Though formally at the level a"^ they are still essential, 
as we know, since they produce the same leading 1/Q effects. This is the Nason-Seymour prob- 
lem of non-inclusiveness of jet observables and it is a part of the whole programme of giving 
unambiguous predictions for the magnitudes of power-behaving contributions. 

Our solution can be formulated as follows. First we give a definite prescription for the 
naive contribution which treats gluon decays inclusively and incorporates the running coupling 
and finite gluon mass effects. To this end we have chosen to substitute \/k^ + for kt both 
in the kinematics and in the definition of the observable, so as to treat a "massive gluon" 
contribution. Given this prescription, we have demonstrated that the two-loop corrections to 
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the naive treatment amount to a universal, that is observable-independent, normahsation effect 
(Milan factor). 

The universality of the Milan factor has three ingredients. Firstly, it relies on universality of 
soft radiation, the latter being responsible for confinement effects. Secondly, it is based on the 
concept of universality of the QCD interaction strength, all the way down to small momentum 
scales, which has been processed through the machinery of the dispersive approach. Finally, it 
includes a certain geometric universality of the observables under consideration. To clarify the 
last point, we note that the contribution of a massless parton z to a given observable V (linear 
in the limit of small parton momenta) can be factorised as 

Vi = hi ■ f^^Hvi) , 

with kti the modulus of transverse momentum and rji the rapidity of the parton, and / a function 
that depends on the observable. For example, for the thrust and C parameter we have 

/(^)(r/) = e-^G(r/) + e^e{-r,) , f^^\r,) = ^ . 

The rapidity integral of this function enters linearly into the magnitude of the power contribution 
and determines the observable-dependent coefficient of the latter. 

The crucial observation is that this very same integral appears both in the terms driven by a 
parent gluon (naive TZq and the inclusive correction 7?.'") as well as in the separate contributions 
from offspring partons the gluon decays into (the non- inclusive correction 7?."*). The two decay 
partons give 

which contribution has to be weighted with the parent gluon decay probability. First we observe 
that the decay matrix element depends only on the relative rapidity Arj = r/i — 772 (Lorcntz 
invariance). Then, that Arj stays finite in the essential region (IR-safety of V). Therefore, 
keeping Ar] and ku fixed, we can perform an overall rapidity integration extending it to infinity: 

Vi + V2 =^ {hi +h2)- p^^^ ■ 

As a result, the dependence on the observable factors out thus ensuring the universality of the 
non-inclusive correction. What is left is the standard integral of the sum of two transverse 
momenta weighted with the two-parton decay matrix element. 

For some specific observables sTich as the height of the transverse momentum distribution 
of Drell-Yan pairs at qt 0, EEC with back-to-back kinematics, etc. we expect an interplay 
between the perturbative resummation and the genuine confinement effects to lead to calculable 
fractional powers of Q in the non-perturbative correction terms. Whether this expectation is cor- 
rect is an intriguing question which should be addressed both theoretically and experimentally. 
The magnitude of these fractional power contributions is determined by the non-perturbative 
parameters Ai and ^2,1 (see (4.1)). The latter parameter enters also in the 1/Q^ power correc- 
tions to DIS (e'^e^) structure (fragmentation) functions, DIS sum rules, the Drell-Yan i^-factor, 
as well as in the mean value of the 3-jet resolution parameter j/a [6,29]. From the analysis of F2, 
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Dasgupta and Webber [29] derived a phenomenological estimate ^2,1 ~ —0.2 GeV^. This implies 
that confinement effects should dampen the height of the Drell-Yan distribution at = 0. 
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